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ABSTRACT 


A Computer program to determine the ultimate failure 
load of a circular bulkhead with horizontal and vertical 
stiffeners iS presented. The circular plate is treated as a 
Grillage with the plating acting as an effective flange on 
the stiffener. The grillage 1S composed of a combination of 
transverse and longitudinal members, however, the spacing 
between nodal points is not uniform. The plate-stiffener 
combination 1S subjected to a uniform hydrostatic pressure 
that is modelled as equal point-forces at each of the 
interior nodes. A matrix system of inequality constraints 
(on maximum moment values) and nodal point equilibrium 
equations is established and the applied point-force is 
maximized using a Linear Programming maximization routine 
(ZX3LP) from the International Mathematical and Statistical 
Library (IMSL). A detailed explanation of the computer 
program 1S presented. 


The ultimate failure load for the unsupported segments 
of plating (between stiffeners) is investigated. A brief 
review of the use of the Upper and Lower Bound Theorems is 
presented. 


Results from the computer program are compared to 
experimental data for similarly stiffened circular plating 
Subjected to a uniform pressure until collapse. 


Thesis Supervisor: Dr. Joao Manuel Gomes de Oliveira 


Title: Associate Professor of Ocean Engineering 
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CHAPTER 1 INTRODUCTION 


BACKGROUND 


In the last decade, the demand for new and more 
efficient modes of marine transportation has forced 
naval architects and marine engineers to search for 
practical and reliable methods of structural analysis 
which would take full advantage of the inherent 
Strength of the materials of construction in providing 


a safe design with a minimal Structural weight. 


Elastic design methods provide an accurate and 
reliable technique et Dredlcting Structural 
deformations resulting from applied design loadings 
whenever the material behavior is linear elastic. 
However, 1n applications where the total deformation is 
not a limiting factor, elastic analysis methods do not 
provide an effective technique for prediction of the 


ultimate failure load of a given structural design. 
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In an attempt to take full advantage of the 
inherent ductility of structural steels, a design 
methodology was developed which treats load levels 
exceeding that necessary to reach initial yield stress 
(elastic limit). This methodology, called plastic 
limit analysis, 1s not useful for predicting local 
Stress levels or deformations, but is a reliable and 
effective method of determining the ultimate collapse 
load. In those cases where structural deformation is 
not a limiting factor, or where the ultimate 
load-bearing capacity of a given structure is a primary 
concern, plastic limit methods of analysis provide the 
Structural engineer with an effective tool for 
accurately determining the load limit of a given 
Structure. If the designer requires additional 
information, such as local stresses and deformations, 
or 1f material Strain hardening or finite deflections 
are an important consideration, then a numerical 
analysis uSing an incremental theory of elasticity will 


need to be performed {1}. 
A detailed review or available literature 
concerning plastic limit analysis of structures with 


stiffened and unstiffened plating indicated a great 
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deal Of recent interest in applying "plastic" 
techniques to Structural analysis and design, with a 
view toward structural optimization. By far, the 
Greatest majority of stiffened-plating applications 
developed utilizing plastic limit analysis dealt with a 
class SE problems that could most generally be 
categorized as stiffened rectangular plating with 


uniformly spaced stiffeners. 


An unstiffened circular plate has been analyzed 
using plastic limit analysis {2} with a linear 
programming technique and the "Lower Bound" Theorem (to 
be discussed in detail). Similarly, a stiffened 
Circular plate was studied using elastic analysis to 
optimize the plate-stiffener combination {3}. In this 
analysis, however, the stiffeners were a radial and 
circumferential system of stiffeners (axisymmetric) 
which were compatible with the axisymmetric circular 
plate, and the stiffened circular plate was optimized 


by maximizing the load per unit deflection. 
In the field of submersible vehicle design, a 
cylindrical hull form provides the best compromise 


between competing criteria of adequate hull strength, 
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maximum internal volume, and minimum hydrodynamic drag. 


As a result, Submersible vehicle bulkheads are 
essentially stiffened circular plates. Due to the 
requirement Oe minimizing fabrication cost, the 


bulkhead stiffening system 1S composed of a grid 
network of horizontal and vertical stiffeners, not 
radial and circumferential stiffeners. Thus, the 
as-constructed submersible vehicle bulkhead cannot be 
analyzed by an adaptation of available elastic (or 
plastic) axisymmetric, orthotropic plate analysis 
methods. Also, ultimate failure is the design 


Criteria, without any limit on deformation. 


In their paperj]{4}, Palermo and Bart applied the 
concept of plastic analysis to the as-constructed 
Submersible vehicle bulkhead. The method detailed in 
this paper combined the concept of a plastic hinge and 
elastic grillage analysis to determine the required 
Scantlings for the main horizontal girder. However, 
the vertical stiffeners were analyzed individually with 
the main horizontal girder assumed to be a fixed 
Support for each half of the vertical stiffener. More 
recent methods allow for plastic analysis of the 


grillage as a whole. 
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This thesis will attempt to analyze the stiffened 
Circular bulkhead as a grillage, adapting plastic 
Grillage limit analysis techniques developed by Hodge 
{5} for a rectangular grillage. The net result will be 
an interactive computer program (CIRCPLAT) designed to 
be inexpensive to run and simple to use, providing a 
reasonably accurate, yet conservative estimate of the 
collapse load of a user-defined, stiffened, circular 


plate. 


The results of the CIRCPLAT analysis will be 


compared to actual bulkhead failure test data. 
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CHAPTER 2 PROBLEM DEFINITION 


FORMULATION OF THE PROBLEM 


The basic problem under consideration is _ to 
determine the failure load of a built-in, circular, 
stiffened bulkhead uSing horizontal and vertical 
Stiffeners when it 1S Subjected to a uniform, lateral, 
hydrostatic pressure. To find an accurate solution in 
the simplest, most efficient manner, it is first 
necessary to decide what is the most effective way to 
model the problem. Reference {6} provided a recent, 
well researched compendium of current methods of gross 
panel analyses along with the inherent assumptions and 
limitations of each method, and a comparison of their 
relative accuracy. Figure 1 provides a brief summary 
of the results of comparison of analyses methods given 
by reference SOs Reference {6} also indicated 
excellent results uSing plastic limit analysis. Based 
on a very favorable comparison between the finite 


element program results and the "grillage analysis", 
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and, the close correlation between both methods and 
experimental data, at appeared to be highly 
advantageous to model the stiffened circular plate as a 
grillage of intersecting beams with an effective width 
of plating acting as a flange to the stiffener. With 
this model, it should be possible to develop a computer 
program that is relatively inexpensive to run, yet 
Sacrifices little or nothing me) accuracy for 


determination of the failure load. 
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Total Stress (psi) Percentage 
Centerline Difference 
Longitudinal 

Stiffener at Fixed *% 


Boundary End 


CrEnoteTop lc ist Order Zo nse 0 Sis 
Plate 
Analysis 2nd Order Zo 505 SatlOe: 
Beam & Plate Ind) Paes) +.04% 
Finite Elements 
Element 
SOLID SAP Orehotropic DBI os) Ve See 
Plate 
Elements 
Finite Element 26 Feo See 
ICES STRUDL 
Grillage Analysis ZS 2. -2% 
Simple Beam Theory 30,184 eal ape: 


** Percentage difference with respect to ICES STRUDL Results 


Bigure 1. 


GRILLAGE MODEL: 


ASSUMPTIONS 


Stiffened Panel Analysis Method Comparison. 


The use of the grillage model in conjunction with 


Chapter 2 


Problem Definition 


14 





plastic limit analysis involves several approximating 


assumptions: 


1. The in-plane compressive loading is small compared 
to the lateral loading and its effect 1s 
negligible. 


2. The influence of the shear force on the formation 
of the plastic hinge is negligible. 


3. The torsional rigidity of the stiffeners is small 
and its effect is negligible. 


4. The effect of the plating can be accounted for by 
incorporating an “effective breadth" of plating as 
a flange on the stiffener. 

5. Plane sections remain plane in the plastic range. 


6. Deflections are infintessimal. 


7. The material is elastic-perfectly plastic; 


8. Drucker's Postulate is valid. 


Reference] {5} discusses the validity of assuming 
that the effect of in-plane stresses are negligible. 
Experimental evidence indicates that if the in-plane 
Stress 1s less than 15 percent of the compressive yield 
Stress, the effect of the axial force can be neglected. 
Current SMCRERISTE pressure hull design, utilizing a 
ring-stiffened cylinder, provides adequate support 


against hull compressive effects sufficient to ensure 
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that compressive axial stresses in the bulkhead 
Stiffeners will not exceed 15 percent of yield stress. 
References {4} and {6} contain correction factors that 
can be applied to the stiffeners' plastic section 
moduli to account for the effect of in-plane stresses 
above 15 percent of yield. These correction factors 
were not incorporated in the CIRCPLAT program. The 
program could easily be modified to include these 
corrections if large axial forces were to be included 


in the analysis. 


The assumption that shear force is negligible is a 
reasonable one for long, slender beams. For short, 
deep beams, a correction factor must be applied to the 
Section modulus. To study the relative magnitude of 
shear and bending stress in a beam, consider a 
rectangular, prismatic beam with a single concentrated 


load as shown in Figure 2: (Reference {6}) 
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Wo 


a 


Figure 2. Rectangular Beam with Concentrated Load. 


The maximum bending moment, which occurs at 
midspan 1S: M=PL/4. The maximum bending stress will be 
at the same point: 


Oimax = Mc/I = (3/2) (PL/bh<). 


The shear force 1S constant at P/2 between each 
Support andthe load P. The formula for maximum shear 
Stress in a rectangular beam is: 

Tmax=3/4)(P/bh) . 
Then Omax/Tmax = 2L/h. Thus, we can see that when L is 


many times (10 or more) greater than h, the magnitude 
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of the maximum shear stress as compared to the maximum 


bending stress 1S negligible. 


Reference]{7} discusses the practical significance 
of torsional rigidity in grillage analysis, and 
concludes that the effect of torsion on the plastic 
moment for I-beams 1s 0.1 percent, and is therefore 


negligible. 


The question of incorporating an "effective" 
breadth of plating as a flange on the stiffener 15S one 
that appears to have a number of answers as to how much 
breadth to include as "effective." Reference {4} used 
a value equal to the stiffener spacing, up to a maximum 
oe S0t? where t=plate thickness. Reference {7} 
utilized a formula developed by Schade: 

Beene eaoy 
where B 1s the fraction of stiffener spacing that is 
effective, and r 1s the ratio of the distance between 
points of zero bending moment and the plate's breadth. 
In reference {8}, extenSive experimental research by 
Kendrick was performed specifically to determine the 
effective breadth of plating to be included as a flange 


to the stiffener when calculating the plastic section 
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modulus. The results of this research showed that 
experimental plastic moments agreed quite closely with 
the theoretic values calculated by assuming all the 
plating effective as a flange to the stiffener. For 
this reason, the plastic section modulus calculation 
incoporated in the CIRCPLAT program utilizes the full 
Stiffener spacing as the "effective" breadth, with no 


maximum limit (such as 60t). 


The remaining assumptions are the standard 
assumptions utilized in deriving the equations for 


plastic limit analysis. 


PLATE MODEL ASSUMPTIONS 


As an additional consideration, it is necessary to 
consider the possibility of plate failure between 
stiffeners. Because this 1s a circular. plate 
intersected by horizontal and vertical stiffeners, the 
resulting sections of unsupported plate created come in 


a variety of shapes (see figure 3). Plastic limit 


Chapter 2 Problem Definition ILS 





analysis can similarly be applied to these plate 


sections to determine their ultimate failure load. 





Baigure 3. Stiffened Circular Plate Geometry. 


In applying plastic limit analysis to thin plates 


it 1S assumed: 
1. Elastic-perfectly plastic material 


2. No in-plane loads 


3. Infinitessimal deflections 
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4, No Shear effects 


5. Johansen's yield surface applies 


Reference {9} details the use of Upper and Lower Bound 
Theorems of Plasticity for determining plate failure 


load subject to the above assumptions. 
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CHAPTER 3 STRUCTURAL ANALYSIS 


BASIC PRINCIPLES 


Reference {1} presents an excellent development of 
the theory of plastic analysis of structures. In 
understanding plastic failure analysis, it is first 
necessary to understand the concept of the plastic 
hinge. Consider pure bending of a prismatic beam of 
rectangular cross-section subject to an applied moment. 
Coordinate axes and sign convention are positive as 


indicated in figure 4. 
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Figure 4. Rectangular Beam subject to pure bending 


The bending moment is aeveu by: 
were f 2Oi(2)9 
-H 
The maximum elastic bending moment will occur when 
the stress at the outer fiber reaches yield. Then the 
moment in this case becomes: 


Me=(4/3) (BH Oo) 


If the applied moment is increased still further, 
the outer fibers will begin to deform plastically, and 


the elastic to plastic interface will begin to move 
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toward the neutral axis until, eventually, k=0 (See 


figure 5.) 
0% ke 
KH 
Figure 5. Elastic-Plastic Stress Distribution 


At the point where k=0, the beam can not support 
any further increase in applied moment, and a plastic 
hinge 1S said to have developed. This fully plastic 
moment, denoted by Mo, for a rectangular section, is: 


Mo = 200BH 


In the partially plastic range between Me and Mo, 


the deformation 1s theoretically controlled by the 


assumption that plane sections remain plane. 
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A yield hinge, then, has the following properties: 


1. If the applied moment is less than Mo, the yield 
hinge will transmit the full moment without 


allowing any rotation. 


2. If the applied moment 1S equal to Mo in magnitude, 
the yield hinge will transmit the full moment, but 
will permit rotations of arbitrary magnitude in the 


direction of the applied moment. 


3. The yield hinge will not transmit any moment 


greater than Mo in magnitude. 


With the concept of a yield hinge defined, failure 
of a structure occurs when sufficient yield hinges are 
formed to allow infinitessimal motion of all or part of 
the structure. The existance of sufficient yield 


hinges 1S termed a "mechanism". 


Several important concepts involved in the idea of 
a failure mechanism are exemplified in the case of a 
Simple, indeterminate structure under a distributed 


load. Consider the case of a beam fixed at one end, 
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Simply Supported at the other end, and under a 


uniformly distributed load P, as in figure 6. 


Oe 2 (a) 


Figure 6. Beam with Distributed Load 


Rivst vig can be seen that the problem is 
indeterminate, because there are 3 unknown reactions 
and only 2 equations of equilibrium (horizontal force 
equilibrium is implicit). From elastic analysis, we 
know that the shape of the moment distribution is that 
of figure 6.c. Thus, we know that a plastic hinge will 
form at the fixed end, and one will form somewhere near 
the center. As soon as the first hinge 1s formed, the 


problem becomes determinate Since we have one less 
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unknown, leaving 2 unknowns and 2 equilibrium 


equations. 


For example; af the first hinge forms at the fixed 
end, then summing moments about point A gives: 
“Mo+R L~PL {2 = 0 
Summing y-directed forces gives: 


RTR, = PL 


These two equations can be solved for both R, and 
R,- Although the problem is now determinate, we do not 
yet have a mechanism. The plastic hinge at A would 
allow rotation, but the simple support at B still 
prevents it, so a second hinge must form somewhere near 
the center of the beam. These two hinges now form a 
collapse mechanism. Plastic hinge rotation can proceed 


unimpeded. 


UPPER BOUND THEOREM 


Continuing with the example of Figure 6, assume 
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that the second hinge forms exactly at the center of 
the beam. The internal work of deformation is the 
product of the plastic moment times the angle through 
which it rotates. For the case of the beam in Figure 
6, if the hinge at the fixed end rotates by 6, the 
center hinge rotates by 28. Work is independent of 
rection cf rotation. Then 


Wi=Mo0+M,(26)=3M,0. 


External work is the product of applied force times 
the distance moved. For a distributed load, the 
average distance moved is 1/4(L°6). If the applied 
load is large enough to do sufficient external work to 
equal the required internal work of the collapse 
mechanism, then collapse will occur. Therefore, 
[P/4\(L* 9 )=3Me will yield the value of P required for 
our assumed mechanism. This result, padieMaye is the 
required distributed load to cause the "assumed" 
collapse mechanism to occur. If the assumed collapse 
mechanism were the exact collapse mechanism, the value 
for P would be the exact value to cause collapse. 
Looking at Figure 6c., the moment distribution shows 
that the maximum moment occurs to the right of the 


center of the beam. As a result, the actual plastic 
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hinge will form to the right of the center. Asa 
result, the collapse load, P, determined above, 1S not 
the exact collapse load. The Upper Bound Theorem 
essentially states that the failure load determined by 
an assumed collapse mechanism that 1S not the exact 
collapse mechanism will be an upper bound to the exact 
failure load. TiuSees ene lOad, plsi2M./L , 1S greater 


than the actual collapse load, and 1s an Upper Bound. 


LOWER BOUND THEOREM 


Remembering that once the first plastic hinge 
forms, the problem becomes determinate, the following 


equations can be written: 


M(Q) = R,L = pie = -Mo 
M(Z) = R,(L-Z)- P(L-2)72 = Mo 
M'(Z) = -Rp + P(L-Z) = 0 


where Z is the point at which the second plastic hinge 
forms. With 3 equations and 3 unknowns, the reaction 
R, distance Z and load P can be determined. The result 


1s: 
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R = 4.8284 Mo/L 
P = 11.65685 M/L“ 
Z = 0.58579 L 


From the moment curve of Figure 6c, it can be seen 
that if the moment at x = 0 1S -Mo and the moment at 


x=0.58579* L 1s +Mo, then the moment can not exceed Mo 


in magnitude anywhere on the _ beam. This, then, 
constitutes a Statically admissible moment 
ese rr buci1on. The Lower Bound Theorem essentially 
States that Le a Statically admissible moment 


distribution can be found that satisfies equilibrium 
and the boundary conditions, and does not violate the 
yield criterion (in this case, M(x)<Mo for all x), then 
the load P that corresponds to this moment distribution 
is a lower bound to the exact solution. For the 
example of Figure 6, the value of P=11.65685 is a 
lower bound, but in this case it also corresponds to 
the exact solution. If the upper bound load, P, were 
determined at the same point on the beam, the same P 
would result, p =P, indicating the exact solution was 


obtained. 
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EFFECTIVE PLASTIC MOMENT OF A STIFFENER 


Several formulations are available to calculate the 
effective plastic moment of a given I-beam, yet the 
problem of incorporating a general formula in a 
computer program for an I-beam, where the user 
Specifies the thickness and breadth of the upper flange 
(plate thickness and stiffener spacing) Separately from 
the remaining Stiffener Scantlings, presents an 
interesting subtlety. The program must. be able to 
distinguish between the case where the neutral exis 
lies in the web and the case where the neutral axis 
lies in the plate. A Simple, accurate, and effective 
technique to accomplish this 1s given in reference {7}. 


Typical stiffened plate geometry is shown in Figure 7. 
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Uy 
Figure 7. Stiffened Plate Geometry for Plastic Section 


Modulus 


When ri tr,>i, the plastic neutral axis lies in the 
web and the formula for the plastic section modulus,Z, 
1S: 


2 
a=B Dt {ri (r ter, +2)ter ry (r,ti)-(r,-1) }/4ar,. 
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When r,+r,<i, the plastic neutral axis lies in the 
upper flange (plate) and the plastic section modulus is 
given by: 

Z 
2=B Dt, {er t4rotratiter, +er,)-r3(ry+ra) Zee / oie 

From chapter two, the value of effective breadth to 
be incorporated in the calculation of plastic section 
modulus was shown to be the full stiffener spacing {8}, 


Without an arbitrary cutoff value. CIRCPLAT sets B 


equal to stiffener spacing for all calculations. 


GRILLAGE COLLAPSE LOAD 


In reference {10}, Hodge applied the Upper Bound 
Theorem to a rectangular grid (hereafter called a 
"grillage") in order to determine its collapse load. 
In the method development, it was assumed that 
concentrated loads were applied at the joints (or 
nodes), of the grillage, although the method could be 
made applicable to distributed loads. Having assumed a 


collapse mode and determined the corresponding upper 
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bound collapse load, it 1S necessary to construct a 
Statically admissible bending moment distribution to 
verify that the assumed collapse mode is feasible. If 
there are no other collapse modes that would yield a 
lower collapse load and the moment distribution 1s 
Statically admissible (the absolute value of all 
moments are at, or lower than, the plastic yield 
moment), then the correct collapse load has been 
determined. Reference {10} contains numerous examples 


of the above procedure. 


This method can be carried out with little 
difficulty for regular rectangular grillages, where the 
loads at all nodes are equal. A regular rectangular 
grillage iS one where the spacing is uniform between 
longitudinal beams, and the spacing is also uniform 
between transverse beams. In cases where these 
requirements are not met, it becomes difficult to 
predict a deformation pattern that will result in an 
improved upper bound. Moreover, the assumed 
deformation pattern may not have a  Sstatically 


admissible moment distribution. 
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Reference {5}, by Hodge and Malone, presented a 
feasible method to overcome these difficulties in 
obtaining a correct collapse load. The proposed method 
involves taking a Lower Bound approach using Linear 
Programming. The load associated with a statically 
admissible moment distribution is, by the Lower Bound 
Theorem, a lower bound to the true collapse load. If 
the load can be maximized while the moment distribution 
1S maintained statically admissible, the maximum lower 
bound load should be the correct collapse load. Linear 
Programming, explained in Appendix 8B, provides a 
Simple, computer method to maximize the load, subject 
to the constraints of static admissibility and the 


equations of equilibrium. 


The grillage of Figure 8 will be discussed in 


detail as an example of the use of the Lower Bound 


Theorem in conjunction with Linear Programming. 
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(a) Grillage with load (b) Symmetric node 


arrangement 


(c) Detail of 


Node B 


Figure 8. 2x2 Grillage with Non-Uniform Loading. 

For this problem the beams are simply supported at 
all boundaries. The plastic yield moment of the 
horizontal and vertical beams 1S Mo. Moments are 
positive when they produce compression in the top 
fibers. 

Chapter 3 Structural Analysis 36 





From Figure 8b we can see that there is symmetry 
about a line through nodes A and C. This Symmetry can 
be used to more effectively solve the problem. As a 
result of symmetry, the number of unknown moments is 
reduced from 8 to 4. The numbers on the ends of the 


beam spans indicate the unknown moments in Figure 8b. 


Because the problem of plastic collapse is only 
concerned with the statically admissible moment field, 
Support reactlons are not required to be solved for 


and, aS a result, are not listed as unknowns. 


Since the loads are applied at the nodes, there are 
no midspan loads and so for each span, the shear forces 
at the ends must be equal and opposite. Moment 
equilibrium at node B for each span means that 
V, =(M;-0)/L, Vo (M4-M), )/L, etc, from Figure 8c. The 
Sum of the shear forces at each node must equal the 
applied load at that node. For node B, the result is: 


(2M.,+2M .-M)-M, )/L = 0 


5 


Or 2M, * 2M 5- MM, = | 


Applying Similar reasoning to the other two 


independent nodes; 
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4M, ~2M,-P = 


\ 
© 


Notice that there are four unknown moments, and 


only three independent equilibrium equations. 


In addition to the equations of equilibrium, the 
linear programming problem requires the equations of 
Somseraint. Each of the moments 1S constrained to be 
equal to or less than the plastic yield moment, 
however, this moment can be positive or negative. The 
linear programming constraint values must all be 
positive. In order to meet this requirement, a new 
variable 1S created, X, =M, +Mo 1=1,2,3,4 so that while 
~Mo<M .<+Mo, OSX .S2Mo. Therefore, Xs 1S constrained to 
positive values. The collapse load is the fifth 
unknown, called X_. 


5 


Substituting M, =X, —Mo into the three equilibrium 


equations the result is: 


eS Els Eee 2, = 2Mo 
4X) ~2%,, ore = 2Mo 
“2K, +k), = 2Mo 
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The load, which is the function to be maximized can be 
defined by a separate equation (called the objective 
function) : 


ns 


By the above procedure, all of the requirements for 
establishing a linear programming problem have been met. 
The next step is to format the information to facilitate 
input to the Linear Programming routine. The required input 
format is a function of the linear program to be utilized. 
At Massachusetts Institute of Technology, one of the 
Standard linear programs available in the International 
Mathematical and Statistical Library (IMSL) 1s called ZX3LP 
(see reference {14}). This program was utilized by CIRCPLAT 


and 1S explained fully in Appendix B. 


ZX3LP requires the input to be arranged in 3 matrices. 
Matrix A contains the coefficients of the unknown variables 
in the equations of constraint, followed by the coefficients 
of the unknown variables in the equilibrium equations. 
Matrix B contains the values from the right-hand side of the 
constraint equations, followed by the values from the right 
hand side of the equilibrium equations. Matrix C contains 


the coefficients of the unknown variables in the objective 
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Bune t. On. Matrices A,B, and C for the sample problem of 


Figure 8 are shown in Figure 9. 


A Matrix B Matrix 


Hl 0 0 0 0 2Mo 
0 1 0 0 0 2Mo 
0 0 1) 0 0 2Mo 
0 0 0 1 0 2Mo 
sedi 2 2 tf 0 2Mo 
4 mz 0 0 = 2Mo 


| 2Mo 


— 
© 
© 

| 
NO 
> 
© 

{_ — 

r- 

bs 


C Matrix [0 0 0 0 1 | 


Figure Se Linear Program ZX3LP Input Matrices. 
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The primary solution of the ZX3LP program is contained 
in array PSOL(N). The values in array PSOL(N) are the 
values of Ky Xo Xz peer eee Sy For the problem of Figure 8, 
the values of PSOL(1) through PSOL(4) must be adjusted by 
the value Mo to give the actual nodal moments. This is 
Beeauscce M; =X. —Mo- The value of the collapse load 1s given 


by PSOL(5) for this case. 


REPLACEMENT THEOREM 


The derivation of a method for determining the collapse 
load of a grillage has, up to now, been for grillages with 
concentrated loads at beam intersections. In extending this 
method to encompass the problem of a uniformly loaded 
Circular plate, it is necessary to discuss the correlation 
between the concentrated load problem and the case of the 
M~corm load. Reference {10} contains a discussion and 
proof of the Replacement Theorem. ‘The essence of the 
Replacement Theorem is that i1f a given beam span of constant 
plastic yield moment, Mo, is acted upon by a distributed 


load whose direction does not change sign on that span, the 
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Span iS not made stronger if the distributed load is 
replaced by equipollent concentrated loads. That is, the 
collapse load of a span under a set of concentrated loads 
will be less than that for the same span under the 


equipollent distributed load. 


An equipollent force system is defined in reference {12} 


as one in which: 


1. The resultant of the forces of one system 1S equal to 


the resultant of forces of the other system. 


2. The sum of the moments of one system about an arbitrary 
base point is equal to the sum of the moments of the 


other system about the same point. 


Using the above theorem, if the uniform load on the 
grillage 1S replaced by an equipollent system of 
concentrated forces at nodal points, then the method of 
grillage analysis just developed, which computes a "best" 
lower bound collapse pressure for concentrated forces, can 
be applied directly to compute a conservative "best" lower 


bound for a uniformly distributed load. 
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Looking at the section of stiffened plate in Figure 10, 
it can be seen that the pressure P distributed over an area 
a x bcan be replaced by a concentrated force of magnitude 


P(a x b). 





baigure 10: Section of a Stiffened Plate 


If the spacing is uniform in the x-direction and 1s also 
uniform (not necessarily the same) in the y-direction, then 
the area a x b will be the product of the stiffener spacing 
in the x andy direction. From this, it can be seen that 
for a regular rectangular Pesce, the concentrated force 
at every node will be the same. In the case of a grillage 


stiffened circular plate, where the stiffeners are equally 
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spaced along the main diameters (x and y directions), this 
Same replacement process works exactly the same for nodes 
not near the circumference of the circle. The replacement 
process for nodes near the circumference is not exact, as 
can be seen in Figure 11. There are nodes where the 
replacement force will be too high, and nodes where the 
replacement force will be too low, but overall, the area of 
the circle 1s greater than the area associated with all of 


the nodes. 





miguge, 14. Top View of a Stiffened Circular Plate. 
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It iS necessary to ensure that the estimate provided by 
the application OF the Replacement Theorem remains 
conservative for the case of the stiffened circular plate. 
In CIRCPLAT, the forces constrained to be equal at all of 
the nodes. Then, after the collapse load is determined, the 
collapse pressure is calculated by dividing the sum of the 
nodal forces by the total area of the circular plate 
[P=(Sum of F)/A . Since the area in the denominator is 
Greater than that associated with the nodal forces (the 
rectangle of Figure i1.), the pressure will be conservative. 
The accuracy of this method is improved as the number of 


beams in the grillage is increased. 


UNSUPPORTED PLATE COLLAPSE LOAD. 


Calculating the collapse load of the grillage does not 
provide any information about the possibility of collapse of 
the sections of plating between stiffeners. To be complete, 
analysis of the stiffened plate must include a determination 
of the collapse pressure of the unSupported plate segments. 


Figure 12 indicates the variety of shapes possible for 
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unsupported segments of a stiffened circular plate. In 
order to simplify the collapse load calculation, use is made 
of a corollary to the Lower Bound Theorem discussed in 
reference {13}. This corollary states that you cannot 
decrease the collapse load of a structure by increasing the 


Strength of any part. 





Figure 12. Unsupported Plate Segments of Stiffened 


Circular Plate. 


In Figure 12 we can see that if the rectangular sections 


of unsupported plating are a x b, the unsupported sections 
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along the perimeter of the circle are less than a x b. 
These sections can be considered to be rectangles of 
dimensions a x b with a wall built into the rectangle for 
increased support. Here, invoking the corollary to the 
Lower Bound Theorem discussed above, it can be seen that 
these rectangles strengthened by walls cannot have a lower 


collapse load than the full rectangular section. 


Reference {1} contains a conServative formula for the 
lower bound collapse load of a clamped rectangular plate 
Subject to a uniform pressure. The unsupported sections of 
plate are best modelled as clamped because the rotation at 


the support 1S zero. 


This formula, implemented in CIRCPLAT, is: 
P = (16M, /b“) (1+8*) 
where a=plate length 
b=plate breadth 
S=b/a 


This equation is a lower bound associated with the 


Johansen yield curve. 
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CHAPTER 4 CIRCPLAT PROGRAM 


PROGRAM DESCRIPTION AND USE 


CIRCPLAT was written using FORTRAN-77 programming 
language on a VAX-i1i computer with the VAX/VMS operating 
system. The program was designed to be run interactively 
with a minimum of required information from the program 
user, while providing the user with the choice of inputting 
a greater amount of detailed information if it 1s considered 


necessary. 


CIRCPLAT was designed to handle a maximum of twenty 
Stiffeners (horizontal plus vertical) but this capacity 
could be increased by modifying the DIMENSION statement at 
the beginning of the program. Similarly, CIRCPLAT can 
accommodate up to a maximum of 12 stiffeners in any one 


direction (horizontal or vertical). 


The stiffener spacing is assumed by CIRCPLAT to be 


uniform along the horizontal and vertical diameters. Thus, 
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when the user inputs the total number of vertical 
stiffeners, the total number of horizontal stiffeners, and 
the bulkhead diameter, the program computes the x and y 
coordinates of every beam intersection, including the 
coordinates of beam-boundary intersections. These 
coordinates are used to determine the span lengths that are 
needed for calculation of the local shear forces in the 


nodal equilibrium equations. 


CIRCPLAT takes advantage of the symmetry that exists 
about the horizontal and vertical diameters and reduces the 
remaining input required (and the calculations) to only 
those concerned with the upper left quadrant of the circle 
(which includes any beam(s) that lie directly on _ the 
horizontal and/or vertical diameters). During user input, 
the horizontal beams are numbered 1,2,3,etc., starting from 
PiemetoOp, counting down to (and including) the horizontal 
beam at the centerline (if there is one). Similarly, 


vertical beams are numbered 1,2,3,etc. from left to right. 


Inputs for the material yield stress, plate thickness, 
then the stiffener scantlings are requested by the program 
and the plastic yield moments are calculated. As a 


convenience, the plastic moment resulting from the first set 


Chapter 4 CIRCPLAT Program 49 





of vertical stiffener scantlings is assigned to all vertical 
Stiffeners, the user can then modify this value, stiffener 
by stiffener, as required. The user has the option of 


revising the plastic moment of each stiffener by either 


1. Inputting a full set of revised scantlings for that 
Stiffener from which the program calculates the plastic 


yield moment; or 


2. Inputting the plastic yield moment value directly for 


the specified beam. 


Once the vertical beams have been fully specified, the 


mentical routine 1S run for the: horizontal beams. 


The option of specifying the plastic yield moment value 
directly for an individual beam was added to allow the user 
to insert a calculated plastic moment corrected for the 
effect of shear stress or axial force if it should become 
necesSary. The method of obtaining these corrections is 


discussed in detail in reference {4}. 


CIRCPLAT was designed for analysis of a grillage where 


all beams are prismatic. For those cases involving beams 


Chapter 4 CIRCPLAT Program 50 





that vary in cross-section along their length, the minimum 
plastic yield moment should be used to obtain a conservative 
failure load, or, the program itself will have to be 
modified to allow for changes in stiffener geometry along a 


given beam length. 


For those beams where shear deformations or effects must 
be taken into account, if the corrected plastic yield moment 
is less than the minimum plastic yield moment for the 
uncorrected beam, the corrected value should be input. 
Otherwise, to remain conservative, the minimum uncorrected 


plastic yield moment should continue to be used. 


The current version of CIRCPLAT is designed to analyze 
the built-in submersible buikhead. The program treats the 


grillage as being fully clamped at all boundaries. 


At this point, the user input is completed, and the program 
begins to calculate the parameters that it will have to 
Paput ee) the ZX3LP linear programming subroutine. 
BaSically, the remainder of the program is designed to 


@ecomplish four functions: 
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1. Loading the correct values into the proper locations in 
Matrices A,B,and C as indicated in the example of 


Chapter Three. 


2. Establishing the exact matrix DIMENSION values required 


by ZX3LP. 


3. Calculating the failure load of the unsupported sections 


of plating between stiffeners. 


4, Formatting the solution produced by ZX3LP to produce a 


properly labelled output. 


The solution is stored automatically ina data file 
labelled FORQ18.DAT. This data file contains a listing of 
the input data and the output data.The output data consists 


of a listing of the following items: 


1. The failure load of the weakest segment of unsupported 


plating. 
2. The failure load of the grillage as a whole. This value 
1s unrelated to the value in (1) above (i.e. they are 


treated as Separate problems). 
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3. The normalized value (with the corresponding plastic 
moment value used to normalize it) of the moment present 
at failure at each of the vertical beam nodes. This is 
followed by a similar listing for all horizontal nodes. 
The listing of node moments is only for those nodes in 
the upper left quadrant of the plate. The rest are 
symmetric. Also, the vertical nodes are numbered from 
the top of the beam downward, working from left to 
right. In a similar fashion, the horizontal nodes are 


listed from left to right, working from top to bottom. 
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CHAP iER 5 ReESUbTS AND DISCUSSION 


COMPARISON OF RESULTS 


The results of CIRCPLAT were compared to the results 
obtained from an actual scale model failure test of a 
stiffened circular bulkhead (Reference {16}). The ratio of 
CIRCPLAT predicted collapse pressure to experimental 
collapse pressure was 0.61. There are a number of valid 
reasons why the CIRCPLAT prediction was 39% lower than the 


experimental value. 


First, the horizontal and vertical stiffeners on the 
bulkhead under test were not prismatic. That is, the 
scantlings varied significantly along the stiffener length. 
To accommodate possibly large shear effects on the main 
horizontal girder, the web depth increased ina linear 
fashion (over approximately one quater span) at each end to 
a web depth almost twice that of the enter Span. 
Variations of web depth of a similiar magnitude were also 


common on vertical stiffeners. Because the program assumes 
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prismatic beams, and to ensure conservative results, the 


smallest value for web depth was input for every stiffener. 


Secondly, the bulkhead design tested incorporated both 
HY-80 and HY-130 steels. Again, because the program allowed 
BOE only one material yield strength, and to ensure 
conservative results only the lower yield stress was used in 


the analysis. 


Finally, a less important factor, the program can 
presently handle a maximum of twelve stiffeners in any one 
direction. There are twelve tertiary horizontal stiffeners 
and one main horizontal girder for the tested bulkhead (for 
a total of thirteen). To maintain the odd-symmetry 
relationship, and keep the total under twelve, it was 
necessary to reduce the number of tertiary stiffeners to 
ten. This has a minor effect, as indicated by the fact that 
the effect of tertiary stiffeners was ignored entirely in 
the analysis of reference {4}. Again, the results will be 


more conservative. 
As a result of the above three factors, CIRCPLAT results 
can be assumed to be very conservative; this is borne out by 


the above comparison. 
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The previous results were for a strengthened version of 
an existing circular bulkhead design. For the original (not 
strengthened) version of this bulkhead, scantling data was 
readily available, but the given collapse load data was 
based on design values, rather than actual failure testing. 
CIRCPLAT was run and results were compared to the design 
collapse load for the circular bulkhead of the original 
design. When CIRCPLAT was run with tertiary stiffeners not 
included, the ratio of CIRCPLAT collapse pressure to design 
collapse pressure was 0.83. When CIRCPLAT was run with all 
tertiary stiffeners (8 of them) included, the ratio was 
eS Js The correlation here is considered good, and the 
effect of tertiary stiffeners for this case would seem to be 
Srqmificant. Failure load analysis of the unsupported 
plating segments for the original design indicated that the 
unsupported plating was the weakest portion. The ratio of 
theoretical collapse pressure to design collapse pressure 
for the unsupported plate (clamped rectangular plate) was 
0.73 for an Upper Bound solution, and 0.49 for a Lower Bound 
se0t10Nn. This indicated the bulkhead plating would not 
withstand full design pressure. This problem did not show 


up in the improved bulkhead design. 
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In reference {10}, Hodge provides a formula for the 
exact collapse load of a uniform circular plate, with a 
fully clamped boundary. The formula reduces to 


Da SER Be fae 


where P uniform collapse pressure 
Oo = material yield stress 
t = plate thickness, 


R = Radius of plate. 


For a plate 10 feet in diameter, 2 inches thick, with a 
yield stress of 50,000 psi, the exact uniform collapse 


pressure 1S 156 psi. 


Referring to the formulas of Chapter 3 for calculating 
the plastic section modulus of a stiffened plate, CIRCPLAT 
can be used to accurately model a circular plate as a 
grillage by setting all web and flange scantlings equal to 
zero, except for the web depth. The web depth can be any 
arbitrary value, as it factors out of the plastic section 
modulus formula, but it must not be zero, because it 1S in 
the denominator of 3 and ry and would cause these values to 
become singular. For convenience, input web depth to 


CIRCPLAT as 1.0, and input all other web and flange 


Scantlings as zero. As a result, the plate will be modelled 
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by CIRCPLAT aS a egrillage of stiffeners with rectangular 
cross section (bxh) where b 1s stiffener Spacing, determined 
by the number of stiffeners the user inputs (more stiffeners 


will give better results); and h is plating thickness. 


EOL the Same plate analyzed above using Hodges's 
formula, CIRCPLAT calculated a failure load of 154.7 psi, 
only 1% less than the exact value. The plate was modelled 
for this case as ten vertical and ten horizontal stiffeners. 
This result 1S significant in that it points out that 
CIRCPLAT does yield a conservative yet very close, lower 
bound solution, in a case where the input values can closely 
approximate the exact physical problem, and also verifies 
that a uniform load can be approximated quite well by point 


loads at all nodes. 


SUGGESTED IMPROVEMENTS TO CIRCPLAT 


In order to reduce the amount of conservatism indicated 
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by the comparison of CIRCPLAT results to actual stiffened 
plate failure test data, several improvements to the program 


are recommended for implementation: 


1. Modify the input routine to provide the user with the 
ability to input a different set of horizontal and/or 
vertical stiffener scantlings at each nodal point. This 
would enable CIRCPLAT to more closely approximate a 


non-prismatic beam; 


2. Modify the input routine to provide the user with the 
ability to input a different yield stress value for each 
stiffener, or for the circular plate. This would, 
again, provide the user with more flexibility in 


modelling a real, physical problem; 


3. Modify the program to include a variable "fixity" 
factor, £, at boundary nodes, ranging in value from zero 
EOP One . This would allow the user to model a simple 
Support (f=0.), a fixed support (f=1.), or some value in 


between. 
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CONCLUSION 


In conclusion, it 1s felt that the linear programming 
method of grillage analysis, based on the Lower Bound 
Theorem, provides a simple, yet powerful tool for accurately 
determining the collapse load of a given structure. 
CIRCPLAT was designed for the problem of the circular plate 
Stiffened by a grillage network, with fixed (clamped) 
boundaries. However, the applicability of this technique 
can be extended very easily to all types of grillage 
collapse analysis. It is strongly felt that this analysis 


tool yields results that are consistent with empirical data. 
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APPENDIX A PROGRAM LISTING 





Noe of Verte Stiffeners 
NeOwoOt Hore stilrreners 
Plate Diameter 

Plate Thickness 
Material Yield Stress 


INPUT 
Sscantlings of most common 


Vertical Stiffener 
(Web Thickness/Depth, 


Flange Thickness/Width) 









Do any 
Vertical 
Stiifeners differ 
From above 
Value ? 


Choose to: 
1. INPUT new 
Scantlings. 
Or 
2- INPUT new 
Plastic Yield 
Mioment. 


















NO 


NEU or a2), (B) 


Repeat process for above 


Horizontal Stiffeners 
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Calculate Failure 
Load of Unsupported 
Plate 









Cali subroutine ZX3LP 
to solve Linear Programming 
Problem 







Store OUTPUT in file 
FOROL8.DAT 





Program Listing 
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10 
ZO 


40 
50 
60 
70 
80 
90 
100 
mae) 
120 
0 
140 
150 
160 
170 
180 
ene, 
200 
21:0 
22.0 
230 
240 
Zoe 
Z.60 
270 
2 80 
Z90 
S18 
BLO 
B20 
5.3 0 
340 
B50 
3 60 
570 
3 80 
B90 
400 
410 
420 
430 
440 
450 
460 
470 
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— — — 


DIMENSION XCOORD(8,8),YCOORD(8,8),VERTMO(8), 
HORZMO(8) ,MV(8,8),MH(8, 8),NNODE(40),F(40), 
NVERT(40) ,MLIN(80) ,MMAX(80) ,C(80),AA(120, 80), 
B(100),HORZL(40), VERTL(40) , KHORZ(40),KVERT(40), 
MPV(45),MPH(45) ,NVLEN(45),IFLAG(45) ,PSOL(102), 
A(120,80),DSOL(120),RW( 8000) ,1W(304) 


REAL HORZMO,MV,MH,MMAX,MLIN,HORZL,MPV ,MPH,HFLAG,MP 


C 

C USS heel RACE RVELYS INPUTS BEAM AND PLATE PARAMETERS FOR 
C Cit eoh erent REGULAR PIEATE OF INTEREST. SYMMETRY ABOUT 
C THE HORIZONTAL AND VERTICAL CENTERLINES IS ASSUMED. 

C INPUT PARAMETERS AND CALCULATIONS, THEN, ARE FOR THE 

C UPPER ober “QUADRANT OF THE CIRCULAR PLATE. 

C 

G254507 


PRtliewem | THIS PROGRAM WAS DESIGNED FOR A LIMIT OF” 
PRG os: PE NEve TOGA “STRERFENERS (HORIZ.+VERT.) 
PRN ie TWELVE MAX. IN ONE DIRECTION.’ 


Y=SKIP(5) 

PRINT *, ’ THIS PROGRAM ASSUMES THAT HORIZONTAL AND’ 

PRINT *, ’ VERTICAL STIFFENERS ARE EQUALLY SPACED ALONG’ 
PRINT *, ’ THE HORIZONTAL AND VERTICAL CENTERLINES, AND’ 
PRINT *, ’ ARE T-STIFFENERS WITH AN EFFECTIVE PLATE WIDTH’ 


>] 
PRINT *, ’ MAKING THEM I-STIFFENERS.’ 
Y=SKIP(5) 
PRINT *, ’ ENTER THE NO. OF VERTICAL STIFFENERS’ 
READ *, N2 
Y=SKIP(5) 
PRINT *, ’ ENTER THE NO. OF HORIZONTAL STIFFENERS’ 
READ *, M2 
Y=SKIP(5) 
PRINT *, ’ ENTER THE DIAMETER (FEET) OF THE CIRC. PLATE’ 
READ *, DIAM 
D=DIAM*12. 
R=D/2. 
WRITE(18, 4000) 


4000 FORMAT(1X,° INPUT INFORMATION’ /) 


C234567 


WRITE(18,4001) N2 


wo l PORMAL@IxX, NO. OF VERTICAL STIFFENERS:’ ,3X,13 /) 


WRITE(18,4002) M2 


HoGeme ORMAT GluX, NO. OF HORIZONTAL STIFFENERS:’ ,3X,13 /) 


WRITE(18, 4003) DIAM 





480 MUG MEORMA Nalwe PLATESDTAMETER? {3X,F6.2,°CFEET)” /) 


490 C 

50) CemcCAlCCUnATE SIME Ae (HORITZ.) “AND Y (VERT.) COORDINATES OF 
Se C NODES (BEAM INTERSECTIONS OR BOUNDARIES) 
52 0 Cc 

53 0 C 

540 YDIST=D/REAL(M2+1 ) 

518 XDIST=D/REAL( N21 ) 

560 C CHECK ODD OR EVEN SYMMETRY 

570 VARI=REALCN2Z )/2.-IFIXC(REAL(N2)/2.)-.1 
580 VAR2 =REAL(M2 )/2 .-LFIX( REAL(M2 )/2.)-.1 
590 LODDV=0 

600 LODDH=0 

610 MECVARL.GT.0)  LrODDV=1 

620 IF(VAR2.GT.0) ILODDH=1 

630 BECLODDV.EO.0) GO TO 50 

640 KCOL=(N2+1)/2+1 

650 VFLAG=0.0 

660 Gor 10 51 

670 50 KCOL=N2 /2+1 

6 80 VFLAG=1.0 

690 oe! DECLODDH.EBOQ.0) GO. TO 55 

700 KROW=(M2+1)/2+1 

710 HFLAG=0.0 

720 CO TO) 56 

730 55 KROW=M2 /2+1 

740 HFLAG=1.0 

750 56 CONTINUE 

760 CG 

770 CoekGOL = NO, OF NODES ON LONGEST HORIZ. BEAM 
780 C KROW = NO. OF NODES ON LONGEST VERT. BEAM 
790 Ce (THis £S ENCLUDING BOUNDARY NODES AND IS FOR 
800 GeetHe UPPER LEFT QUADRANT ONLY) 

810 C H/VFLAG = 0 IF ODD SYMMETRY 

820 C = 1 IF EVEN SYMMETRY 

830 C 

840 C234 

850 Poel twrkt Ze ALE COORD. MATRICES TO ZERO 
860 C 

870 DO 70 I=1,KROW 

8 80 DOmoOs J=i KCOL 

890 RCOORD CL ,J)=0.0 

900 MeEOORD CI, J)=0.0 

910 60 CONTINUE 

920 70 CONTINUE 

930 C 

940 Cee LLU erNe COORD. MATRICES 

950 Cc 





960 
970 
980 
B90 
1000 
1010 
O20 
1030 
1040 
1050 
1060 
1070 
1080 
O90 
1100 
el O 
2 0 
1130 
1140 
£150 
P60 
el 7 0 
1180 
9:0 
i200 
i 10 
2 0 
B23 0 
1240 
i250 
1260 
i 70 
2 80 
2 90 
1300 
i 10 
ih 2 0 
Es 3 0 
1340 
1350 
1360 
i370 
13 80 
3 90 
1400 
1410 
1420 
1430 


> 


DO 110 NROW=1,KROW 

Om OOM NIC Ol = Colmer 
PP CNuROWw.Gie1) GOsmOo 115 
XVAL=R-(REAL(KCOL-NCOL)+VFLAG/2.)*XDIST 
YVAL=SQRT(R*R-(R-XVAL) *#2 ) 
YMIN=(REAL(KROW)-2+HFLAG/2.)*YDIST 
TEST1=YVAL-YMIN 
eK uicen OD Gon nO iG 
XCOORD(NROW, NCOL)=XVAL 
YCOORD(NROW,NCOL)=YVAL 
COTO 100 
XMAX=SQRT(R*R-((REAL(KROW-NROW)+HFLAG/2.)*YDIST)*#2 ) 
YMAX=SQRT(R*R-( (REAL(KCOL-NCOL)+VFLAG/2 .)*XDIST) *#2 ) 


C234567 


C 
C 


C 


106 


OS 


108 


1G 


100 
110 


PE CNCOR EO eI mCO TO” f 16 
XVAL=R-(REAL(KCOL-NCOL)+VFLAG/2.)*XDIST 
XCOORD(NROW, NCOL)=XVAL 
TEST2=XMAX-(R-XVAL) 
PCUeGH? SLEOn GO TO) 105 
YVAL=(REAL(KROW-NROW)+HFLAG/2.)*YDIST 
TEST3=YMAX-YVAL 
MuGhMors hn YDEST) CO. TO 110 
PP ChEGts sch 0) Cor TO. 108 

YCOORD(NROW, NCOL)=YVAL 

60 EO 100 

XCOORD(NROW, NCOL)=R-XMAX 

GO TO 106 

YCOORD(NROW, NCOL)=YMAX 

60 20-110 

XCOORD(NROW, NCOL)=R-XMAX 

YCOORD(NROW, NCOL)=(REAL(KROW-NROW)+HFLAG/2.)*YDIST 

CO Wa ie 

CONTINUE 
CONTINUE 


USER INPUTS STIFFENER SCANTLINGS TO DETERMINE PLASTIC MOMENTS 


PRINT *, ’ 
BREN (tee 


INPUT MATERIAL YIELD STRESS (PSI) AND PLATE’ 
THICKNESS (INCHES) AS FOLLOWS: ’ 


PRINT *, ’ SIGMA-Y, THICKNESS’ 
READ *, SIGMAY,T2 
Y=SKIP(5) 


Pie SUE INITIAL SCANTLINGCS YOU INPUT WILL BE ASSIGNED’ 
Poo ee tOPALE VERTICAL STIFFENERS. YOU WILL THEN HAVE THE 
PRINT *, “© OPTION OF REVISING SPECIFIC STIFFENERS BY NUMBER 
Ee Neb (BY ENTERING EITHER SCANTLINGS OR PLASTIC MOMENT) ’ 
NCSI ILE (Coy 

IVERT=0 
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1440 
1450 
1460 
1470 
1480 
1490 
1500 
n> 10 
mp2 0 
1530 
1540 
i > 0 
1560 
B70 
> 60 
9 0 
1600 
6 10 
1620 
63 0 
1640 
1650 
1660 
1670 
1680 
ina 9'0 
7 00 
7 10 
1720 
H730 
1740 
my 50 
1760 
in? / 0 
1780 
ih 90 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
ine 7 0 
1880 
ie 90 
1900 
e910 


1, 


4006 FORMAT (1X, 
ie WE Rati CK NES Gin 
11X,’WEB 
ie oeOb UO MC inNGH Sh LOKNESS +“ 
11X,’ BOTTOM FLANGE 


PRINT *, “ INPUT INITIAL VERTICAL BEAM SCANTLINGS AS FOLLOWS: 
PRINT *, “ ENTER WEB THICKNESS,WEB DEPTH, BOTTOM FLANGE ’ 
PRINT *, ’ THICKNESS,BOTTOM FLANGE WIDTH (ALL IN INCHES)’ 

READ *, 73 ,DW,T1,Bl 

LF(LVERT.EQ.0) GO TO 4005 
WRITE(18,4006) NBEAM-1,T3,DW,T1,B1 
‘VERTICAL BEAM NO.:’ 3X,13,/ 
Ses ele INCHES” ./ 
ee ee  INGHES * | / 
,3X,F8.3, 
Beir 8 3. 


DEPTH 
INCHES’, / 
WIDTH INCHES’ ,//) 


GO TO 16 


4005 WRITE(18, US) SEGMAY IZ 13", DW El, B81 


4004 FORMAT(1X, 
Pie MATERIAL YIELD STRESS: 


VERTICAL BEAMS a, 
,3X,F12.2,° PSI’,/ 


11X,°PLATE THICKNESS: ’ ,3X, F8.3, INCHES’, / 
11X,’COMMON WEB THICKNESS: Bye WSs SENCHES  / 

11X,’COMMON WEB DEPTH oes. ENCHES” | / 

le COMMON BOTTOM FLANGE THICKNESS. ot rey INCHES’, / 
11X,’COMMON BOTTOM FLANGE WIDTH SN Boe ENGHES.<20/ 15) 


16 LF(IVERT.EQ-1) GO TO 21 


a0 


DORA OUSI-2-KCOL 
VERTMO(L)= MEGDW XDA S Tere «Ble Th sO LGMAY,T3 ) 
CONTINUE 
COMO 2 
VERTMO(NBEAM)=MP(DW,XDIST whcob ls. 1S LGMAY »T3)) 
PRINT *, “© ARE THERE ANY VERT. BEAMS WHICH ARE DIEEERENT- 
PRON? 7 ‘ FROM THOSE INPUT ABOVE ?’ 


PRINT *, “ OQ=NO, 1=YES’ 
READ *, IVERT 

PEGiVERT. hO.0) GO To 30 
SSK IED) 


PRINT *, ’ 
PRINT *, ’ 
SE 2) 
PRINT *, “© COUNTING IN SEQUENCE LEFT TO RIGHT, INPUT THE’ 
PRINT *, “ BEAM NO. FOR WHICH THERE ARE NEW SCANTLINGS’ 
READ *, NBEAM 

NBEAM=NBEAM+1 

PRINT *, “ DO YOU WANT TO INPUT SCANTLINGS, OR PLASTIC MOMENT: 
PRINT *, ’ O=INPUT PLASTIC MOMENT,1=INPUT SCANTLINGS’ 
READ *, NOPT 

IF(NOPT.EQ.1) GO TO 15 


SINCE SYMMETRY IS ASSUMED, INPUT ONLY THOSE VER 
DEMS SON OR lO THE LEFT OF VERT. CENTERLINE. “ 


PRINT*, “ INPUT PLASTIC MOMENT’ 
READ *, VERTMO(NBEAM) 

N=SKUP(2 ) 

GOmron 22 

LHORZ=0 
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2 0 
1930 
1940 
19 5 0 
1960 
980 
no 9 0 
2000 
2010 
ZO2 0 
Z03 0 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
110 
220 
Zl 3 0 
2140 
eel 5 0 
2160 
eal / 0 
e780 
Za 9 0 
2200 
ee 10 
pe 2 0) 
22350 
2240 
Eee 50 
2260 
2270 
22 80 
E290 
2300 
25 10 
Zo 20 
E30 
2340 
po 50 
E26 0 
23/0 
25 80 
ES 90 
2400 


PRINT *, * NOW A SIMILAR ROUTINE WILL BE FOLLOWED FOR’ 
Peet “aes” HORIZONTAL STIFFENERS.’ 
sal Ne et a INPUT HORIZ. BEAM SCANTLINGS AS FOLLOWS: “ 

3 PRN The ENTER WEB PUEGRNESS.) WEE DEPTH, BOTTOM FLANGE’ 
PRIN 25 <% THICKNESS, BOTTOM FLANGE WIDTH (ALL IN INCHES)’ 


READ *, T3,DW,T1,B1 
IF(IHORZ.EQ.0) GO TO 4007 
WRITE(18,4008) MBEAM-1,T3,DW,T1,B1 

4008 FORMAT(1X,’HORIZONTAL BEAM NiOng no Me 130 
11X, “WEB THICKNESS:’ ,3X,F8.3,’ INCHES’, / 
eer a ee 3k F8.3,° INCHES’ , / 
11X, "BOTTOM FLANGE THICKNESS:’ ,3X,F8.3,’ INCHES’, / 
L1X," BOTTOM FLANGE DEPTH :’,3X,F8.3,’ INCHES’ ,//) 
couro 34 

4007 WRITE(18,4009) T3,DwW,T1,B1 

4009 FORMAT(1X,’ HORIZONTAL BEAMS’, / 
11X, "COMMON WEB THICKNESS: ’ ,3X,F8.3,° INCHES’ ,/ 
11X,°COMMON WEB DEPTH :’,3X,F8.3,’ INCHES’ , / 
11X, COMMON BOTT.FLG.THICKNESS:’ ,3X,F8.3,’ INCHES’ ,/ 
11X, COMMON BOTT.FLG. WIDTH :’,3X,F8.3,’ INCHES <0) 

34 IF(IHORZ.EQ.1) GO TO 3] 
DOs =2 |. KROW 
HORZMO(1I)=MP(DW,YDIST,T2 ,B1,T1,SILGMAY,T3) 

53 CONTINUE 
GO Towd? 

Bal ee EAM =MPE (DW YDIST,T2 ,B1,T1,SI1GMAY,T3) 

32 PRINT *, “ ARE THERE ANY HORIZ. BEAMS WHICH ARE ’ 


PRINT *, “ DIFFERENT FROM THOSE INPUT ABOVE 2° 
PRINT *, ’ O=NO, 1=YES’ 
READ *, IHORZ 
IF( IHORZ.EQ.0) GO TO 40 
= Ser) 

C234567 
PRINT *, “ SINCE SYMMETRY ASSUMED, INPUT ONLY THOSE’ 
PRINT *, “  HORIZ. BEAMS ON OR ABOVE HORIZ. CENTERLINE’ 
Mei LPC?) 
PRINT *, “ COUNTING IN SEQUENCE FROM TOP TO BOTTOM, ’ 


PRINT *, “INPUT THE BEAM NO. FOR THE NEW SCANTLINGS. ’ 
READ *, MBEAM 

MBEAM=MBEAM+1 

PRINT *, * DO YOU WANT TO INPUT SCANTLINGS OR PLASTIC MOMENT?’ 
PRINT *, °  O=INPUT PLASTIC MOMENT, 1=INPUT SCANTLINGS’ 

READ *, MOPT 

GECMOPEsEQ.1) GO TO 35 

PRINT *, “ INPUT PLASTIC MOMENT’ 

READ *, HORZMO(MBEAM) 

Y=SKIP(2) 
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2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
Be 10 
zed 2 O 
> 0 
2540 
Ee 5 0 
2560 
Be / 0 
2580 
Zo 90 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
Bo 7 0 
2680 
2690 
2700 
2710 
2720 
273 0 
2740 
27 50 
2760 
27 70 
80 
a 90 
2800 
2810 
2820 
E63 0 
2840 
2350 
2860 
2870 
2530 


GO TO 32 
40 CONTINUE 


é 
é 
C INITIALIZE ALL MV AND MH FLELDS TO ZERO 
C 
¢ 
DO 90 I=1,KROW 
DO 80 J=1,KCOL 
MV(L,J)=0.0 
MH(1I,J)=0.0 
80 CONTINUE 
90 CONTINUE 
6234567 
C 
C 
C INITIALIZE NODAL EQUATION COUNTER TO ZERO AND COUNT 
C INTERNAL NODES FOR EACH LINE. SUM INTERNAL NODES TO 
C COUNT NODAL EQUATIONS 
C 
C 


NODEEQ=0 
DO 120 J=2,KROW 
XMAX=SQRT(R*R-((REAL(KROW-J)+HFLAG/2 .)*YDIST)*#2 ) 
NNODE(J)=LFIX(XMAX/XDIST-VFLAG/2 .+0.9999) 
NODEEQ=NODEEQ+NNODE(J) 
120 CONTINUE 


C 
C 
C CREATE A ONE DIMENSIONAL ARRAY OF UNIT NODAL GRID LOADS 
C 
e 
DO 130 K=1,NODEEQ 
F(K)=1 
130 CONTINUE 
KKK=1 


ICKMKsel 0) 1BC=2 
WGK Kact. 0) CO TO’.131 
PRINT *, ’ ENTER BOUNDARY CONDITIONS: ’ 
PRINT *, ’ 1=SIMPLE SUPPORT, 2=FIXED SUPPORT’ 
READ *, IBC 
131 CONTINUE 


END SOESUSER  LNPUT 


CC e) 


WRITE(18,4014) 
4014 FORMAT(1X,° OUTPUT INFORMATION ia) 
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2890 
2210 0 
eo 1 0 
p72 0 
ZS 0 
2940 
po 50 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
Be 50 
3060 
B07 0 
3080 
5090 
e100 
3110 
B12 0 
B30 
3140 
B50 
3160 
B70 
B80 
BL 90 
Bz 00 
Bz 10 
B22 0 
Be 30 
3240 
B2 20 
3260 
62 70 
52 80 
52 90 
5300 
Be 10 
a> 2 0 
33 0 
3340 
BS 50 
3360 


PM CNOL tahoe OAND MORTsRO.1sAND.1T3 .LT.0.00001) Go To 4020 
BB=YDIST/2. 

L=XDIST/2. 

iP CxpiGhaGEsYDIST) = BB=xXDIST/2 . 

[IPGOESteCE. YDIST) L=YDIST/2. 
PL=SIGMAY*(T2 **2 /BB**2 )*(1+BB**2 /L**2 ) 

WRITE(18,4015) PL 


Sole tORMA LG en UNaUr es PLATE SECTION COLLAPSE LOAD: “,F10.4) 
C 

C 

C ASSIGN MOMENT FLAGS TO HORIZ. BOUNDARY NODES IAW INPUT B.C. 
C 

C 

C 

€753456/7 


4020 DO 140 NN=2,KROW 


MCOL=KCOL-NNODE(NN) 
MH(NN,MCOL)=1. 
LF(LBC.EQ.1) MH(NN,MCOL)=0.0 


140 CONTINUE 
C 
C COUNT INTERNAL NODES ON VERTICAL BEAMS 
é 
DOM MS Omb—2) KOOL 
YMAX=SQRT(R*R-((REAL( KCOL-L)+VFLAG/2 .)*XDIST) ¥**2 ) 
NVERT(L)=LFIX( YMAX/YDIST-HFLAG/2 .+0.9999) 
150 CONTINUE 
234567 
C 
C ASSIGN MOMENT FLAGS TO VERT. BOUNDARY NODES IAW INPUT B.C. 
C 
C 
DO 160 JK=2 ,KCOL 
MROW=KROW-NVERT (JK) 
MV(MROW,JK)=1. 
LF(IBC.EQ.1) MV(MROW,JK)=0.0 
160 CONTINUE 
© 
C SET INTERNAL NODE MOMENT FLAGS 
6 
DO 180 JL=2 ,KROW 
NSTART=KCOL-NNODE(JL)+1 
DO 170 JN=NSTART,KCOL 
MV(JL,JN)=1. 
MH(JL,JN)=l. 
170 CONTINUE 
180 CONTINUE 
C234567 
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3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 
3710 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3 830 
3840 


C 
C CREATE A LINE ARRAY OF NON-ZERO VERTICAL MOMENTS AND AN 
C ASSOCIATED ARRAY OF MAXIMUM PLASTIC MOMENTS; COUNT 
C TOP TO BOTTOM, LEFT TO RIGHT 
C 
NLV=1 
DO 210 J=2,KCOL 
DO 200 I=1,KROW 
RCM Ipc 60 TO 1200 
MLIN(NLV)=MV(I,J) 
MMAX(NLV)=VERTMO(J) 
NLV=NLV41 
200 CONTINUE 


Zor CONTINUE 


C 

Cee CONTINUE THE LINE ARRAY WITH NON=ZERO HORIZ. MOMENTS AND THEIR 
CemasSOCTATED LIMITS. 
C 

C 


234567 
NLH=NLV 
DO 230 1=2,KROW 
DO 220 J=1,KCOL 
MiCiuGl T)sbr.1l) CO TO 220 
MLIN(NLH)=MH(I,J) 
MMAX(NLH)=HORZMO(L) 
NLH=NLH+1 
220 CONTINUE 
230 CONTINUE 


CALCULATE TOTAL NUMBER OF EQUATIONS (NODAL EQUILIBRIUM PLUS 


INEQUALITY EQUATIONS) 
NEQTOT=NODEEQ+NLH-1 


CoO C 2 


Poole G MATRIX (OBJECTIVE FUNCTION VECTOR) 
Plow lowtne FINAL C MATRIX 


Gr @)"™ C2 


DO 240 I=1,NLH 
Car i=00 
240 CONTINUE 
C(NLH)=1.0 


INITIALLIZE THE AA MATRIX (VARIABLE COEFFICIENTS ) 
NOT eNO ss OF ROWS REQUIRED BY ZX3LP IS NEQTOT+2 


De i a Sal ee ee 
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Bo) 0 
3 860 
3 870 
3 880 
3890 
B70 0 
5710 
3920 
B25 0 
3940 
3950 
3960 
By /0 
3980 
B90 
4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
42 80 
4290 
4300 
4310 
4320 


DO 250 I=l1,NEQTOT+2 
DO 245 J=1,NLH 
AA(I,J)=0.0 


245 CONTINUE 

Zo0 CONTINUE 

#2354567 

C 

ee hal Gea wdt  BeMATR IX (RIGHT HAND SIDE VECTOR) 

C FILL IN ANY KNOWN (CINEQUALITY) VALUES FOR B 

Cee te eSOUNDARIES ARE FIXED, THIS 15 THE FINAL B MATRIX 
C 


DO 255 I=1,NEQTOT+2 
ME@mecEeNER—1) CO TO 254 
B(1)=2.*MMAX(TL) 


CG) ONOY 705) 5 
254 B(1I)=0.0 
Zo) CONTINUE 


C 
C FILL IN AA MATRIX INEQUALITY CONSTRAINT VARIABLE COEFFICIENTS 
C 
DO 265 I=1,NLH-1 
AA(I,I1I)=1. 
265 CONTINUE 
C 
C LOAD NODAL FORCE COEFFICIENTS INTO THE AA MATRIX 
C 
INDEX=NLH 
DO 270 I=1,NODEEQ 
AAC INDEX,NLH)= -1.0 
INDEX=INDEX+1 
270 CONTINUE 


CREATE A LINE ARRAY OF HORIZ. LENGTHS, AND ANOTHER ARRAY 
OF VERTICAL LENGTHS 


HORIZONTAL ARRAY: 


Cl@) Gave) Gy ©) 


NNN=1 
DO 280 J=2 ,KROW 
NSTART=KCOL-NNODE(J)+1 
DO 275 L=NSTART,KCOL 
HORZL(NNN)=XCOORD(J,L)-XCOORD(J,L-1) 
NNN=NNN+1 
275 CONTINUE 
280 CONTINUE 
C 
C VERTICAL ARRAY: 
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4330 
4340 
4350 
4360 
4370 
43 80 
43.90 
4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 


NNN=1 
DO 290 J=2 ,KCOL 
NSTART=KROW-NVERT(J)+1 
DO 285 I=NSTART,KROW 
VERTL(NNN)=YCOORD(I-1,J)-YCOORD(I,J) 
NNN=NNN+41 
285 CONTINUE 
290 CONTINUE 


CREATE AN ARRAY (KHORZ(N)//KVERT(N)) WHICH ASSOCIATES THE 
SEQUENTIAL (HORIZONTAL//VERTICAL) NODE NUMBER TO EACH INTERNAL 
NODE NUMBER. INTERNAL NODE NUMBERING Lit haetcur  ,TOP-BOTTOM 


Se ies fa ae aw ep 1 


NINT=1 
DO 300 I=2 ,KROW 
NSTART=KCOL-NNODE(1L)+1 
DO 295 J=NSTART,KCOL 
IF(IBC.EQ.1) KHORZ(NINT)=NINT 
IF(IBC.EQ.2) KHORZ(NINT)=NINT+(I-1) 
C2345678901234567890123456789012345678901234567890123 4567 8901234567890 
IF(J.EQ.2) KOUNT=0O 
HEV Gle2)) GO TO 2.96 


298 KVERT (NINT) =KOUNT+1I+1IBC-1-(KROW-NVERT(J) ) 
GO TO 299 
296 KOUNT=0 


DO 297 K=2,J-1 
KOUNT=KOUNT+NVERT(K)+IBC-1 


017 CONTINUE 
GO TO 298 
299 NINT=NINT+1 
295 CONTINUE 
300 CONTINUE 
C 
C CREATE A LINEAR ARRAY OF HORIZ. AND VERT. PLASTIC MOMENTS 
C ASSOCIATED WITH EACH INTERNAL NODE 
C 
NINT=1 
DO 310 I=2,KROW 
NSTART=KCOL-NNODE(I)+1 
DO 305 J=NSTART,KCOL 
MPV(NINT)=VERTMO(J) 
MPH(NINT)=HORZMO(LI) 
NINT=NINTH1 
305 CONTINUE 
310 CONTINUE 
C 


C CREATE AN ARRAY ASSOCIATING VERTICAL LENGTH INDICES WITH 
C INTERNAL NODES 
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4810 C 


4820 NINT=1 
4830 DO 330 I=2,KROW 

4840 NSTART=KCOL-NNODE(1)+1 

4850 DOM3 25.) =NSTART, KCO1, 

4860 IF(J.EQ.2) KOUNT=0 

4870 inseseGr a? ) EO) TO: 32 6 

4880 Bee NVLEN(NINT)=KOUNT+I-(KROW-NVERT(J)) 

4890 GOe lo 325 

4900 326 KOUNT=0 

4910 DO 327 Ran. 7Si 

4920 KOUNT=KOUNT+NVERT(K) 

4930 B27 CONTINUE 

4940 GO) oe 6 

4950 329 NINT=NINTH1 

4960 Be5 CONTINUE 

4970 330 CONTINUE 

mgs0 60. 

4990 C CREATE AN ARRAY OF FLAGS ASSOCIATED WITH EACH INTERNAL NODE 
5000 C TO INDICATE THE FOLLOWING FLAG CONDITION: 

Oe (1) NODE IS NEAR VERTICAL BOUNDARY 

O20 C (2) NODE IS NEAR HORIZONTAL BOUNDARY 

5030 (3) NODE LS NEAR BOTH BOUNDARIES 

5040 C (4) NODE IS A SYMMETRY NODE ABOUT THE HORIZ. DIAMETER 

5050 Cc (5) NODE IS A SYMMETRY NODE ABOUT THE VERT. DIAMETER 

S060 C (6) NODE IS NEAR BOUNDARY AND NEAR HORIZONTAL AXIS OF SYMMETRY 
mo70 C (7) NODE IS NEAR BOUNDARY AND NEAR VERTICAL AXIS OF SYMMETRY 
BOS0 C 

5090 NINT=1 

5100 DO 320 I=2,KROW 

1 0 NSTART=KCOL-NNODE(1)+1 

B20 DO 315 J=NSTART,KCOL 

5130 IF(HORZL(NINT).LT.XDIST) ILFLAG(NINT)=2 

5140 IF(VERTL(NVLEN(NINT)).LT.YDIST) IFLAG(NINT)=1 

5150 IF (VERTL(NVLEN(NINT)).LT.YDIST.AND.HORZL(NINT).LT.XDIST) 
2160 = =€234567890123 4567890123 4567890123 4567890123 456789012345678901234567890 
5170 1 IFLAG(NINT)=3 

5180 IF(J.EQ.KCOL) IFLAG(NINT)=5 

5190 iC tOnkeOlwANDensnOn2) LRUACCNINT)=7 

5200 IF(L.EQ.KROW) IFLAG(NINT)=4 

5210 IF(I.EQ.KROW.AND.J.EQ.NSTART) IFLAG(NINT)=6 

5220 NINT=NINT+H1 

5230 315 CONTINUE 

5240 320 CONTINUE 

5250 DO 1000 I=l1,NODEEQ 

5260 WRITE(6,1001) HORZL(1),NVLEN(1),VERTL(NVLEN(L)) 

5270 1000 CONTINUE 

52 80 1001 FORMAT(1X,F8.3,18,F8.3) 
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5290 
53 00 
5310 
5320 
5330 
5340 
5350 
53 60 
B3 /0 
53 80 
Bo 90 
9400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
9480 
5490 
5500 
B10 
a2 0 
BD 3 0 
2540 
pb» 5 0 
E> 00 
2 gy 
Do SO 
BS 29 0 
3600 
5610 
poz 0 
5630 
5640 
9650 
5660 
5670 
56 80 
3690 
5700 
by 10 
9/20 
5730 
5740 
5) oye 
5760 


C234567 


C 


C LOAD NODAL EQUATION COEFFICIENTS INTO AA MATRIX 


C 


343 


345 


INDEX=NLH-1 
DO 340 I=1,NODEEQ 


IF(IBC.EQ.1) GO TO 341 
KHVAL=NLV-1+KHORZ(L) 

KVVAL=KVERT(1) 

IF(I.EQ.NODEEQ) GO TO 342 

AAC INDEX+1,KHVAL-1)= -1./HORZL(1) 
IF(IFLAG(1).EQ.5.AND.IODDV.EQ.1) GO TO 343 
PEGUNEAGEL)- EO. 5, AND .LODDV.E0.0) GO TO 344 
IF(IBC.EQ.2.AND.IFLAG(1).EQ.7.AND.IODDV.EQ.1) GO TO 343 
IF(IBC.EQ.2.AND.IFLAG(I).EQ.7.AND.IODDV.EQ.0) GO TO 344 
AACINDEX+I,KHVAL)= 1./HORZL(1I)+ 1./HORZL(I+1) 

AAC INDEX+I1,KHVAL+1)= -1./HORZL(I+1) 

COemoO 345 

AAC INDEX+I,KHVAL-1)= -2./HORZL(1) 

AAC INDEX+1,KHVAL)= 2./HORZL(1) 

Gor To s45 

AAC INDEX+I,KHVAL-1)= -1./HORZL(L) 

AACINDEX+1,KHVAL)= 1./HORZL(TI) 

KVVAL=KVERT(L) 

IF(L.EQ.NODEEQ) GO TO 342 

AACINDEX+1,KVVAL-1)= -1./VERTL(NVLEN(L)) 

SCE CACGE).EO.4]4ND, LODDH.EO.1) GO To 346 

Se GLQE cl) -f054 .AND.TODDH.EO.0) €O TO 347 
IF(IBC.EQ.2.AND.IFLAG(1).EQ.6.AND.IODDH.EQ.1) GO TO 346 
IF(IBC.EQ.2.AND.IFLAG(I).EQ.6.AND.IODDH.EQ.0) GO TO 347 
AACINDEX+I,KVVAL)= 1./VERTL(NVLEN(L))+ 1./VERTL(NVLEN(1)41) 


C23 45678901234567890123 45678901234567890123456789012345678901234567890 


346 


347 


342 


348 


349 


AACINDEX+1,KVVAL+1)= -1./VERTL(NVLEN(L)41) 
GO TO 340 

AACINDEX+1,KVVAL-1)= -2./VERTL(NVLEN(L)) 
AACINDEX+1,KVVAL)= 2./VERTL(NVLEN(L)) 

GO TO 340 

AACINDEX+I,KVVAL-1)= ~1./VERTL(NVLEN(L)) 
AACINDEX+1I,KVVAL)= 1./VERTL(NVLEN(L)) 

GO TO 340 

IF(IODDV.EQ.1)GO TO 348 

AAC INDEX+I,KHVAL)= 1./HORZL(1I) 

AAC INDEX+I,KHVAL-1)= -1./HORZL(L) 

GO TO 349 

AACINDEX+I,KHVAL)= 2./HORZL(L) 

AAC INDEX+1,KHVAL-1)= -2./HORZL(1I) 
IF(IODDH.EQ.0) GO TO 350 

AAC INDEX+1,KVVAL)= 1./VERTL(NVLEN(L)) 


Te 





57/0 
5780 
Be 90 
5800 
po 10 
Bro 2 0 
5830 
5840 
5850 
5860 
5870 
5880 
Bro 9 0 
5900 
Be 10 
B20 
5930 
5940 
pe 0 
5960 
5970 
5980 
90 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
6110 
oi2 0 
6130 
6140 
6150 
6160 
5170 
6180 
6190 
6200 
a2 10 
6220 
6230 
6240 


(De al ep Jil Shas Ge Bak @ 


S000 
21010: 35 
341 


AAC INDEX+I,KVVAL-1)= -1./VERTL(NVLEN(L)) 
G0n £01340 
AA(INDEX+1,KVVAL)= 2./VERTL(NVLEN(I)) 
AAC INDEX+1,KVVAL-1)= -2./VERTL(NVLEN(I)) 
CONTINUE 
PRCRGeE Osco) TOS 4 I 
WRITE(18,3000) 
Roe CON eC AACI J) d=1,NLH).L=1 , NEQTOT+2 ) 
FORMAT(1X,’THE AA MATRIX IS:’) 
PORMAL CIN 13°02 Ky 44) ) 
PRINT *, ’ SIMPLE SUPPORT NOT FULLY IMPLEMENTED YET’ 
PVCrBe.bO.1) STOP 
M1=NLH-1 
M2 =NODEEQ 
N=NLH 
TA=M2+M1+2 
LRW1=IA*IA+3 *M1+2 *M24+4 
IW1=2 *M2+3 *M1+4 
MN=MAXO(N,M1+M2 ) 
CALL DIMEN(A,AA,B,C,N,M1,M2,1A,IRW1,IW1,PSOL,DSOL,IW,RW,MN) 


Cony Rie 2.sbe SOLUTION VECTOR PSOL(I) TO THE CORRESPONDING 


NORMALIZED MOMENTS AT THE NODES, 


AND ALSO DETERMINE THE 


APPLIED NODAL FORCE AND EQUIVALENT UNIFORM PRESSURE 


400 


500 


600 
450 


DO 400 I=1,NLH-1 
MLIN(I)=(PSOL(I)-MMAX(1I))/MMAX(IL) 
CONTINUE 
PECGHPLAGsEOQ. 1.) 
LF(HFLAG.EQ.0) 
LF(VFLAG.EQ.1) 
UP CV RDAG. EO .0) 
FNODE=PSOL(NLH) 
FTOT=FNODE*(4.*NODEEQ-HFAC#2 .*NNODE(KROW)-IFAC*2 .*NVERT(KCOL) 


HFAC=0 
HFAC=1 
IFAC=0 
bien Ca 


ietHEAG*LEAC) 


PRESS=FTOT/(3.14159265*R¥R) 
WRITE(18,450) PRESS 
WRITE(18, 460) 
WRITE(18,470) 
DO 500 I=1,NLV-1 
WRITE(18,480) I,MLIN(L),MMAX(TI) 
CONTINUE 
WRITE(18, 490) 
WRITE(18,485) 
DO 600 I=NLV,NLH-1 
WRITE(18,480) I-NLV+1,MLIN(I) ,MMAX(I) 
CONTINUE 
FORMAT(1X,’THE COLLAPSE PRESSURE IS:’,F10.2,’ PSI’) 
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460 FORMAT(1X,’THE NORMALIZED VERTICAL MOMENTS ARE: 
1(TOP-BOTTOM, LEFT-RIGHT)’ ) 
470 FORMAT(1X,’VERT MOMENT NO.’,5X,’NORMALIZED VALUE’ ,5X, 
1’MAX PLASTIC MOMENT’ ) 
480 FORMAT(8X,13,9X,F17.4,5X,E15.4) 
485 FORMAT(1X,°HORZ MOMENT NO.’,5X,’NORMALIZED VALUE’ ,5X, 
1’MAX PLASTIC MOMENT’ ) 
490 FORMAT(1X,’THE NORMALIZED HORIZONTAL MOMENTS ARE: 
1 (LEFT-RIGHT, TOP-BOTTOM)’ ) 
6234567890123 4567890123 4567890123 456789012345678901234567890123456789( 


510.2 
END 


C234567 


SUBROUTINE DIMEN(A,AA,B,C,N,M1L,M2,1A,IRW1,IW1,PSOL, 
1 DSOL,IW,RW,MN) 
DIMENSION AA(120,80) ,A(LA,N),B(IA) ,C(N),PSOL(MN) 
DIMENSION DSOL(IA),IW(IW1),RW(IRW1) 
DO 12 J=1,N 
DO 11 I=1,1A-2 
A(I,J)=AA(I,J) 
iL CONTINUE 
12 CONTINUE 
CALL ZX3LP(A,1IA,B,C,N,M1,M2,S,PSOL,DSOL,RW, IW, IER) 
RETURN 
END 


THIS SUBPROGRAM CALCULATES PLASTIC MOMENT OF EQUIVALENT 
SPEPPEaNER. WHEN CALLING THIS FUNCTION. SPACE” IS THE 
STIFFENER SPACING BETWEEN ADJACENT STIFFENERS. 


(Die Se a Sp a | 


Ge 





100 C 


110 FUNCTION MP(DW,SPACE,T2,Bl1,T1,SIGMAY,T3) 

2 0 REAL MP 

0 Rl CoWwst3)/ CSPACHE*12 ) 

140 R2=( B= El) / CSPACE*T2 ) 

150 R3=T2 /DW 

160 R4=T1/DW 

170 LFCRI+R2.LT.0.99999) GO TO 5 

180 Z=SPACE*DW*T2*(R1*¥(R1+2.*¥R34+2.)4+2.*RIL*R2* (R441. )-(R2-1.)**2 ) 
190 1/(4.*R1) 

200 GO sO: 6 

210 5 Z=SPACE*DW*T2 *(2 .*R1+4.*R2+R3*(1.42.*R14+2.*R2)-R3 *(RI+R2)**2 
220 1+2.*R2*R4)/4. 

230 6 MP=Z*SIGMAY 

240 RETURN 

250 END 

260 c 

270 C 

2 80 C THIS SUBPROGRAM CREATES BLANK LINES ON SCREEN 
290 C 

300 c 

210 

320 O234567 

330 FUNCTION SKIP(N) 

340 DOet0e t—1,N 

350 SKIP=I1 

360 PRINT *, ’ ®’ 


ES 








10 CONTINUE 
RETURN 
END 
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APPENDIZ A 


INTENTIONALLY BLANK 


Program Listing 
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APPENDIX A 


INTENTIONALLY BLANK 


Program Listing 
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APPENDIX B Zio SEN Sg 


ZX3LP LINEAR PROGRAM. 


This section, taken from reference {15}, 1s intended to 


provide a brief explanation of the capabilities and uses of 


the linear programming subroutine ZX3LP. 


IMSL ROUTINE NAME ~- 
PURPOSE = 


USAGE = 


ARGUMENTS A = 


IA = 


APPENDIX B ZX3LP 


7 Seine 

SOLVE THE LINEAR PROGRAMMING PROBLEM 
VIA THE REVISED SIMPLEX ALGORITHM, - 
EASY TO USE VERSION 

CAL) 2xone VAelAP Bae N Mi .M2.S,PSOL, 
DSOL,RW, IW, IER) 

MATRIX OF DIMENSION M1+M2+2 BY N 
CONTAINING THE COEFFICIENTS OF THE Mil 
INEQUALITY CONSTRAINTS IN THE FIRST Mil 
ROWS FOLLOWED BY THE COEFFICIENTS OF 
THE M2 EQUALITY CONSTRAINTS. (INPUT) 
THE LAST TWO ROWS OF A ARE USED ONLY 
AS WORKING STORAGE. 

ROW DIMENSION OF MATRIX A EXACTLY AS 
SPECIFIED IN THE DIMENSION STATEMENT 
IN THE CALLING PROGRAM. (INPUT) 

TWO ROWS OF A ARE REQUIRED FOR WORKING 
STORAGE, AND THEREFORE, IA MUST 
NOT BE LESS THAN M1i+M2+2. 

VECTOR OF LENGTH M1+M2+2 CONTAINING THE 
RIGHT HAND SIDES OF THE INEQUALITY 
CONSTRAINTS IN ITS FIRST M1 LOCATIONS 
FOLLOWED BY THE M2 RIGHT HAND SIDES OF 
THE EQUALITY CONSTRAINTS. (INPUT) 

THE LAST TWO ELEMENTS OF B ARE USED 
AS WORKING STORAGE. 

VECTOR OF LENGTH N CONTAINING THE 

COEFFICIENTS OF THE OBJECTIVE FUNCTION. 
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(INPUT) 


N - NUMBER OF UNKNOWNS IN THE MODEL. (INPUT) 
M1 - NUMBER OF INEQUALITY CONSTRAINTS. (INPUT) 
M2 - NUMBER OF EQUALITY CONSTRAINTS. (INPUT) 
S ae ALU ORs THES OBIECTIVE FUNCTION. (OUTPUT) 


PSOL - VECTOR OF LENGTH N CONTAINING THE PRIMAL 
SOLUTION. (OUTPUT) PSOL IS ALSO USED 
AS WORK STORAGE AND THEREFORE MUST 
HAVE LENGTH AT LEAST MAX(N,M1+M2). 
DSOL - VECTOR OF LENGTH M1+M2+2 CONTAINING THE 
DUAL SOLUTION. (OUTPUT) 
THAT IS, DSOL(1)....,DSOL(M1+M2) 
CONTAIN THE SOLUTION TO THE 
PROBLEM MIN BT*Y SUBJECT TO AT*Y IS 
GREATER THAN OR EQUAL TO C AND Y 
GREATER THAN OR EQUAL TO 0 WHERE 
AT=A-TRANSPOSE AND BT=B-TRANSPOSE. 
WHEN THE PRIMAL PROBLEM HAS 
EQUALITY CONSTRAINTS, THE 
CORRESPONDING COMPONENTS OF THE 
DUAL SOLUTION ARE UNCONSTRAINED. 
DSOL(M1+M2+1) AND DSOL(M1+M2+2) 
ARE USED AS WORKING STORAGE. 


RW - WORK VECTOR OF LENGTH 
(M1+M2+2)*(M1+M2+2)+ 3*M1+2*M2+4., 
IW - WORK VECTOR OF LENGTH 2*M2+3*M1+4. 
IER - ERROR INDICATOR. (OUTPUT) 


TERMINAL ERROR 

IER = 130 INDICATES THAT IA IS LESS 
THAN M1+M2+2. 

IER = 131 INDICATES THAT THE COST 
CRITERION HAS UNBOUNDED VALUES. 

IER = 132 INDICATES THAT THE MAXIMUM 
NUMBER OF ITERATIONS WAS REACHED IN 
ZXOLP. 

IER = 133 INDICATES THAT NO FEASIBLE 
SOLUTION EXISTS. 

WARNING (WITH FIX) 

IER = 70 INDICATES THAT SOME 
ARTIFICIAL VARIABLES REMAINED IN THE 
SOLUTION BASIS AT A ZERO LEVEL AFTER 
PHASE 1. THIS CONDITION CAN BE CAUSED 
BY HAVING REDUNDANT CONSTRAINTS. 
NEVERTHELESS, A SOLUTION IS COMPUTED 
AND RETURNED IN PSOL AND DSOL. 

PRECISION/HARDWARE - SINGLE AND DOUBLE/H32 
- SINGLE/H36,H48,H60 
REQD. IMSL ROUTINES- UERTST,UGETIO, ZXOLP 
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NOTATION - INFORMATION ON SPECIAL NOTATION AND 
CONVENTIONS IS AVAILABLE IN THE MANUAL 
INTRODUCTION OR THROUGH IMSL ROUTINE 
UAE DE 
REMARKS ZX3LP IS INTENDED TO VERY EASY TO USE. THEREFORE. 
THE CALLING PARAMETERS ARE VERY SIMPLE AND THERE ARE 
NO OPTIONS. ZXOLP CAN BE USED DIRECTLY IN SITUATIONS 
THAT REQUIRE MORE FLEXIBILITY THAN IS PROVIDED BY 
Php SAS 


ALGORITHM 


To solve the linear programming problem, 
Maximize c,PSOL, + ...+ CyPSOL, = SP 
Subject to 


aq PSOL + =... + ain PSOL Sb, Pade Ma 
Sepa oOliet  ... 2 ail PoOLy=b 5 ibe Ge lbp »M 
PSOL,=0 iets us sce 
where M = M1 + M>. 
The dual linear programming problem is, 
minimize b,DSOL, + ... + by,DSOLy = SD 
Subject to 
a,,DSOL, + ... + ay 4 DSOLy2c- = leper enn 


DSOL;=0 Dee eis Gy Ml 9 
DSOL; unrestricted in sign when aM ye cea, 


ZX3LP computes the solution to the primal problem, PSOL, 
the solution to the dual problem, DSOL, and the values of 
the objective function S=SP=SD. 


ZX3LP calls ZXOLP which solves the linear programming 
problem by the revised simplex method. 


See reference {14} for additional detail. 


PROGRAMMING NOTES 


ZX3LP is intended to very easy to use. Therefore, the 
calling parameters are very simple and there are no options. 
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ZXOLP can be used directly in situations that require more 


flexibility than is provided by ZxX3LP. 


EXAMPLE 
Suppose we want to Maximize x, + 3x = S&S 
subject to x, = 1 
Xo eS ll 
x} + Xo = el 5 
Fe Nie alee, pa, 
x,20,x5 2 0 
We can proceed as follows: 
INTEGER IA,N,M1,M2,IW(16),IER 
REAL ECGme nb (oO) GZ), RWi 52) -PSOL(4) ,DSOL(6),S 
N = 2 
Mi = 4 
M2 = 0 
IA = 6 
alr 0. 
A = Ge ae 
liv I. 
eas or 
es 
B = ler 
ces 
- 5 
C = Cie oe.) 


CALL ZX3LP (A,IA,B,C,N,M1,M2,S,PSOL,DSOL,RW,IW, IER) 


ty a re 
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APPENDIX C 


SAMPLE OUTPUT 


Sample Output 
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INPUT INFORMATION 

NO. OF VERTICAL STIFFENERS: 10 
NO. OF HORIZONTAL STIFFENERS: 10 
PLATE DIAMETER: 10.00( FEET) 


VERTICAL BEAMS 


MATERIAL YIELD STRESS: JOCC0R I) “P51 

POATLE THICKNESS: 2.000 —INCHES 

COMMON WEB THICKNESS: 020007 INCHES 

COMMON WEB DEPTH : Pe VO0s INCHES 

COMMON BOTTOM FLANGE THICKNESS: 0.000 INCHES 
COMMON BOTTOM FLANGE WIDTH : 02000" INCHES 


HORIZONTAL BEAMS 


COMMON WEB THICKNESS: C2000 PINCHES 
COMMON WEB DEEL Te: Pau OP TNC HES 
COMMON BOTT.FLG.THICKNESS: 0.000 INCHES 
COMMON BOTT.FLG. WIDTH : 0.000 INCHES 


OUTPUT INFORMATION 
Viele Ate oECTION COLLAPSE LOAD: 14/722.2227 


Tne COLLAPSE PRESSURE 18: 4. Oe Bo | 
PteVORMAbE ZED VERTICAL MOMENTS ARE: (TOP-BOTTOM, LEFT-RIGHT) 
Ton MOMENT NO. NORMALIZED VALUE MAX PLASTIC MOMENT 
| =. 0006 0O.35455E+06 
2 O02 0.54 555406 
3 1.0000 0.5455E+06 
4, NE OS ONG, OR 4555 G00 
5 = OOO @ 0.2 bo ONG 
6 -0.0847 OM 2495 E400 
i 1.0000 0.5455E+06 
8 1.0000 0.5455E+06 
, 1.0000 0.5455E+06 
10 = 100 0 0.5455E+06 
Ie) = OUMG 0.5495 E406 
dz 1.0000 024555406 
13 1.0000 0.5455E+06 
14 026/95 0.3455E+06 
gS 1.0000 02 435 E406 
16 =1. 00:00 OPS 455) 5 oe ORS) 
ys =0'..935.6 0. 5455E+06 
18 =. 2.0000 0.5455E+06 
19 -0. 15:66 2545 5E +06 
20 iP. 0008 Hi. 4 55 8+06 
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Zl =O 26024 0.5455E+06 


wD -1.0000 0.5455E+06 
23 -0.3511 0.5455E+06 
24 -0.0847 0.5455E+06 
aE O08 71 0.5455E+06 
26 -0.1988 0.5455E+06 
27 1.0000 0.5455E+06 
THE NORMALIZED HORIZONTAL MOMENTS ARE: (LEFT-RIGHT,TOP-BOTTOM) 
HORZ MOMENT NO. NORMALIZED VALUE MAX PLASTIC MOMENT 

l -1.0000 0.5455E+06 
2 eye 0.5455E+06 
3 1.0000 0.5455E+06 
4 1.0000 0.5455E+06 
5 -1.0000 0.5455E+06 
6 26 283 10 0.5455E+06 
i -0.8532 0.5455E+06 
8 O72 70 0.5455E+06 
9 1.0000 0.5455E+06 
10 -1.0000 0.5455E+06 
11 ; -1.0000 0.5455E+06 
12 -0.6990 0.5455E+06 
13 0.2891 0.5455E+06 
14 1.0000 0.5455E+06 
15 1.0000 0.5455E+06 
16 -1.0000 0.5455E+06 
17 -1.0000 0.5455E+06 
18 -0.4728 0.5455E+06 
19 -0.3433 0.5455E+06 
20 -0.8523 0.5455E+06 
Dl 1.0000 0.5455E+06 
2D -1.0000 0.5455E+06 
23 0.0480 0.5455E+06 
24 Oa7226 0.5455E+06 
25 1.0000 0.5455E+06 
26 1.0000 0.5455E+06 
Da -1.0000 0.5455E+06 





INPUT INFORMATION 

NO.) OF VERTICAL STIFFENERS: 4 
NOw OR WHORTZONTAL STIFFENERS: 5 
PLATE DIAMETER: 27. OO CH Er) 


VERTICAL BEAMS 
MATERIAL YIELD STRESS: 
PLATE THICKNESS: 
COMMON WEB THICKNESS: 
COMMON WEB DEPTH : 
COMMON BOTTOM FLANGE THICKNESS: 
COMMON BOTTOM FLANGE WIDTH 


1.000 INCHES 


HORIZONTAL BEAMS 

COMMON WEB THICKNESS: 
COMMON WEB DEPTH 

COMMON BOTT.FLG.THICKNESS: 
COMMON BOTT.FLG. WIDTH 


OUTPUT INFORMATION 


DNSUPPe er DALE SECTION COLLAPSE LOAD: 
Pie mMooObhArboh PRESSURE IS; a2 243) PSL 
THESNORMALIZED VERTICAL MOMENTS ARE: (TOP~BOTTOM, LEFT-RIGHT) 
MAX PLASTIC MOMENT 


VERT MOMENT NO. NORMALIZED VALUE 
= ee OOO 

=. COO 0 

-0.4645 

1.0000 

=120000 

0.1584 

1.0000 


Oro 02 


OnAnUW & WM 


Tie eORMALIZED HORIZONTAL MOMENTS ARE: (LEFT=RIGHT ,TOP-BOTTOM) 
MAX PLASTIC MOMENT 


HORZ MOMENT NO. NORMALIZED VALUE 
-1.0000 

0.0455 

1.0000 

=. 0000 

1.0000 

ly CO00 

= 15 WOON 

-0.8478 

1.0000 


OMAN A ULF WM He 


SvO0US00 PSI 


02500 INCHES 
5.000 INCHES 
he O00) INCHES 
4.000 INCHES 


Cr OO SEN CHES 

6.000 INCHES 
1.000 INCHES 
4.000 INCHES 


588.3488 


OP Hh Shee cha eI) 
02213 2-807 
O32 32 Bao 
On 2 a ee 
O.2 13 2ETOy 
O72 Va tO, 
022432840 7 
0.21328 207 


0.2450E+07 
0.2450E+0/7 
0.2450E+07 
0.2450E+0/7 
O0.24505+07 
Oy 2 4508207 
O. 24505707 
0.2450E+07 
0.2450E+07 


90 
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